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Abstract. For one-dimensional growth processes we consider the distribution of the 
height above a given point of the substrate and study its scale invariance in the limit 
of large times. We argue that for self-similar growth from a single seed the universal 
distribution is the Tracy- Widom distribution from the theory of random matrices and 
C^ , that for growth from a flat substrate it is some other, only numerically determined 

S ' distribution. In particular, for the polynuclear growth model in the droplet geometry 

the height maps onto the longest increasing subsequence of a random permutation, 
from which the height distribution is identified as the Tracy- Widom distribution. 



PACS numbers: 02.50.Ey, 68.35.Ct 

(N 

E^ . 1. Introduction 

(N 

O I Quite analogous to equilibrium systems close to their critical point, growth processes 

Q^ I exhibit statistical self-similarity and scaling. Their theoretical investigation is a 

continuing enterprise for a variety of reasons. Perhaps the most fundamental one is 

^ ' that statistical self-similarity emerges from the dynamical rules without any particular 

fine-tuning, which is believed to be a basic characteristic of many nonequilibrium 

'^ ' systems |jl|, H, ^. The persisting theoretical challenge is to develop renormalization 

Q ' group type methods which have been so powerful in the understanding of equilibrium 

J-^ ' critical phenomena. 

In the context of KPZ theory Q] most of the recent efforts have concentrated 
/\ ' on scaling exponents and in particular on the question of an upper critical dimension 

^J rl— 1 I— H 

•^ ' B H) 05 m ■ III this note we will focus on the full height distribution at one point above 

the substrate. In addition we will restrict ourselves to one-dimensional growth processes, 
where static fluctuations are the same as in equilibrium due to a fluctuation-dissipation 
property, i.e. the roughness exponent is known to be \. 

The most commonly adopted set-up is to start from a fiat substrate, i.e., if h{x,t) 
denotes the height above x G M at time t, then initially h{x, 0) = ux with slope u. 
h{x,t), t > 0, is random through the dynamical rules which will have to be specified 
for a concrete model. The slope u is a locally conserved quantity. By translation 
invariance, h{x, t) — ux has a distribution independent of x and we may consider the 
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height H(t) = h{0, t) above the origin. For large times it will increase linearly with time 
with a slope dependent growth velocity v{u). H{t) fluctuates on the scale t^/^. Therefore 

H{t) = v{u)t + C{u)t^''^i (1) 

for large t. The random variable ^ is of order one and independent of t and u. Equation 
(|^) should be viewed in analogy to the central limit theorem for a random walk with 
independent increments. In that case v{u) would be the average drift, the fluctuations 
diffusive on the scale t^^^, and ^ a standard Gaussian random variable. 

Following Krug et a/ @, the prefactor C{u) depends only on macroscopic properties 
of the growth model, namely the leading nonlinear part of the slope dependent growth 
velocity, \{u) = v"{u), and the static susceptibility or roughness amplitude, 

A{u) = lim lim x'\{h{x, t) - h{0, t) - uxf), (2) 

X—>00 t^rOO 

with limits in this order. On dimensional grounds we must have then 

C{u)=sign{X{u)){a\\{u)\A{ufY^\ (3) 

where a is a dimensionless number, independent of u, fixing the absolute scale of ^. 
Krug et al ^ have studied ^ on the level of moments for a variety of growth models, 
thereby strongly supporting the claim that the scaling form (|1]) together with (^ is the 
same for all growth models in the KPZ universality class. 

In the very early studies on growth models, like the Eden model, one started with 
a single seed at the origin to have a cluster growing through irreversible aggregation. 
For large times the droplet grows then self-similarly with a shape depending on the 
particular dynamical rule. In fact the macroscopic shape is determined through the 
slope dependent growth velocity by the Wulff construction, i.e. growth velocity and 
shape are related as surface tension and crystal shape in equilibrium. Thus, as for 
a flat substrate, in the droplet geometry there is self-similar growth, but now with a 
(macroscopically) non-vanishing curvature. The analogue of (|l|) is the distance of the 
surface from the origin along a given ray and it is natural to expect for it the scaling 
form 

R{t) =vt + Ct^/^x (4) 

with some universal random variable x- 

The main part of our note studies the shape fluctuations (Q) in the droplet geometry. 
In Section ^ we show that for the polynuclear growth (PNG) model (see for example 
]TT| and references therein) x has the Tracy- Widom (TW) distribution appearing in 
the theory of random matrices. We give some numerical and analytical support for 
the universality hypothesis. In addition, we clearly demonstrate that ^ ^ X- Thus for 
self-similar growth with zero noise initial conditions there are two distinct classes, flat 
substrate (zero curvature) and droplet (non-zero curvature). In both geometries the 
height has the scaling form (|l]), resp. (^, but the universal random variables ^ and x 
have different distributions. 
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2. The PNG droplet 

The PNG model describes a crystal growing layer by layer on a one-dimensional 
substrate through random deposition of particles. They nucleate on existing plateaus 
of the crystal forming new islands. In an idealization these islands spread laterally with 
constant speed by steady condensation of further material at the edges of the islands. 
Adjacent islands of the same level coalesce upon meeting and on top of the new levels 
further islands emerge. 

Nucleation events occur independently and uniformly in space-time. Given the 
space-time coordinates (xo,to), xq G M, to > 0, of such an event, the corresponding 
island nucleates in level 

ho = h{xo,to) + l (5) 

in units of transverse lattice spacings. We also choose the time unit such that the lateral 
growth speed equals one. To determine the height h{x, t) we note that it depends only 
on the nucleation events in the backward light cone of {x,t). We label them as (x„,t„), 
n = 1, 2, . . .. Then 

hn = h{xn,tn) + l, n = 1, 2, . . . , (6) 

according to (^), and h{x,t) is obtained as the maximum over all hn in the backward 
light cone, 

h{x, t) = max{/i„ : n = 1, 2, . . . , |x — x^l < t — t„}, (7) 

with the rule that h{x,t) = if there is no such nucleation event. @ and (0) together 
define recursively all /i„ and therefore h{x,t). 

For droplet growth on the initially flat substrate, a single island starts spreading 
from the origin and further nucleations take place only above this ground layer which 
we define to have zero height. We want to study the probability distribution of h{x,t), 
i.e. Prob{/i(x, t) = n} for n = 0,1,2,..., |x| < t. Clearly h{x,t) is determined 
by the set of nucleation events inside the rectangle R(x,t) = {{x',t') G M^ : \x'\ < 
t' and |a; — x'l < t — t'}. In light-like coordinates, u = {f + x')/\^, f = (t' — a;')/\/2, 
the rectangle R^^^t) equals [0, f/] x [0,1/] with U = {t + x)/V2, V = {t - x)/V2. 
The nucleation events in R(x,t) are Poisson distributed with density one. We label 
them as (M„,f„), n = 1,...,N, N random, such that < Ui < ■ ■ ■ < U]\r < U. 
The corresponding order in the second coordinate v, < fp(i) < ■ ■ ■ < fp(Ar) < V, 
defines then (with probability one) a permutation p of length A^, compare with Figure 
1. In order to distinguish nucleation events corresponding to different height levels, 
we partition {p{l) , . . . , p{N)) into decreasing subsequences according to the following 
algorithm. For the first subsequence one starts with the first entry in the tuple and, 
scanning from left to right, one adds an entry to the subsequence if it is smaller than 
the so far last entry in this subsequence. Every used entry is marked as deleted. When 
the end of the tuple is reached the first subsequence is completed. This procedure 
is repeated to obtain further subsequences by using only the undeleted entries of the 
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original tuple. As a result the permutation p is partitioned into decreasing subsequences. 
The nucleation events of each such subsequence are on the same height level and distinct 
subsequences correspond to distinct heights. The space-time height lines, across which 
h{x',t') increases by one unit, are thus constructed by connecting the nucleation events 
belonging to the same subsequence by a zigzag line, as depicted in Figure 1. In this 




Figure 1. Space-time picture of a tiny PNG droplet 



example we have the permutation (4, 7, 5, 2, 8, 1, 3, 6). The decreasing subsequences are 
(4, 2, 1), (7, 5, 3), (8, 6). They define three height lines. For an arbitrary permutation p, 
h{x,t) equals the number of decreasing subsequences, which is precisely the length, /, 
of the longest increasing subsequence of p. To see this take any increasing subsequence 
of p. By construction every element of this subsequence belongs to a different height 
line. Therefore / < h{x,t). The reverse inequality is obtained by noticing that a 
nucleation event {i,p{i)) at level A; + 1 is necessarily situated in the forward light cone 
of at least one nucleation event (j,p(j)) at level k, i.e. i > j and p{i) > p{j). {j,p{j)) 
is called a predecessor of {i,p{i)). Since there is at least one nucleation event in level 
h{x,t), successive selections of predecessors in consecutive levels result in an increasing 
subsequence of length h{x,t), therefore / > h{x,t). 



In fact, the mapping described here has been used before in [0, where the authors 
consider the Hammersley process []13[ in order to study the asymptotics oi I as N ^ oo. 
We remark that the world lines of particles in the Hammersley process are exactly 
the space-time step lines of the PNG model, after a rotation of 45° together with an 
appropriate adjustment of the boundary conditions. To our knowledge this simple 
observation has not yet been reported in the literature and establishes a connection 
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between the two models in much the same way as between the simple exclusion and the 



zero-range processes |T^ 



For given N the nucleation events are uniformly and independently distributed in 
the rectangle R(x,t) ■ Since the u- and ^-integrations factorize, this induces the uniform 
distribution on the set of all permutations of length A^. By the Poisson law, N has the 
probability distribution 

PiN) = exp(-A) — , X = UV={t'~ x')/2, iV = 0, 1, . . .. (8) 

Thus the height of the PNG droplet is equal to the length of the longest increasing 
subsequence of a random permutation with Poisson distributed length, which for a 



fixed length of the permutation is known as Ulam's problem |]15 . 

In recent years there has been spectacular advance on Ulam's problem, for a 
survey see [1^. One important element are deep combinatorial identities expressing 



the probability distribution of / through quantities known from the theory of random 
matrices. In our particular context the identity reads 

Prob{Z(A) =n} = ^^^^ f exp (2v^X]cos^,) H \^''' " ^^''=1'^"^' (9) 

^•l ^J J[-7r,7r]" j^^ l<j<k<n 

which corresponds to the partition function for the eigenvalues in the unitary ensemble 



with an external cosine potential. Baik et al |T^ succeeded to extract from (^ the 
A ^ oo asymptotics with the results 

/(A)/v^ ^ 2 in probability (10) 

and 

/(A)-2VA .,..,. 

—j^ ^ X 111 distribution, (11) 

together with the convergence of moments. Here x is a random variable, whose 
distribution is the famous Tracy- Widom (TW) distribution of the properly shifted 
and rescaled largest eigenvalue of a GUE distributed random matrix [J18l. x has the 
distribution function 

FW^exp(-f ,.-*).(.)=..), (12) 

where u{x) is the unique solution of the Painleve II equation, 

u" = 2u^ + xu, (13) 

with the asymptotic boundary conditions u{x) ~ Ai(x) as x — >■ oo, Ai(x) being the Airy 



function [jT9l. The probability density, F'{t), of .^ has the asymptotics InF'(t) ^ — |t|^/12 
for t — > — oo and InF'(t) ~ — |t^^^ for t — > cxd. 

We remark that the connection of the right-hand side of equation @ to the Painleve 
II equation has been established on an heuristic level already in the context of unitary- 



matrix models [^ . 
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To come back to the PNG droplet we only have to use that h{x,t) = /(A) and to 
translate (|I0]) and (|TT]) into the surface picture. Equation (|I0|) states that the droplet 
acquires for large times a deterministic ellipsoidal shape given by 

lim t-^h{ct, t) = V2VI - c2. (14) 

i— »oo 

for — 1 < c < 1. Thus with our particular initial conditions the PNG model grows in 
the form of a droplet. The second result ensures that the shape fluctuations are on the 
scale t^/^ and in addition identifies the limiting distribution, 

,. hict,t)-tV2VT^^ ,,^, 

lim = = Y. 15 



Note that in ([15|) only the asymptotic average is subtracted. Therefore x may have a 
mean different from zero, which is indeed the case, as can be seen in Figure 2, where 
the TW distribution F'{t) is plotted on a semi- logarithmic scale. 

3. Flat substrate 

For fiat initial conditions, by translational invariance, it suffices to study the distribution 
of the height above the origin, i.e. H{t) = h{0,t). At first glance, since the droplet 
curvature vanishes on the microscopic scale, one would expect the distributions for ^ in 
Equation (|I]) and for x in Equation @ to be identical. To better understand this point 
we again consider the PNG model. 

The fiat initial conditions in the PNG model are h{x, 0) = and nucleation is 
now allowed everywhere on the line. We consider H{t) = h{0,t). H{t) depends 
on nucleation events in the backward light cone of (0,t), i.e. in the triangle Tt = 
{{x',t'); < t' < t, \x'\ < t}. As before we construct a random permutation from 
the nucleation events inside Tt and H{t) equals the length of the longest increasing 
subsequence. Unfortunately the induced measure on the permutations is no longer 
uniform and the analogue of (|^) does not seem to be known. 

To relate H{t) to the distribution h{0, t) of the height process from the PNG droplet 
we define Hx{t)^ \x\ < t as the height which results if the height lines are constructed 
from nucleation events inside the rectangle {(x',t') : \x'\ <t — f and |a; — x'| < t'} only. 
Since a longest increasing subsequence of nucleation events for H{t) must be contained 
in at least one of these rectangles, we have H{t) = m.a.x\j.\<tHx{t). Now for a given 
set of nucleation events in Tt the mapping t' \-^ t — t', x' ^—>- —x' yields a realization 
of nucleation events for the PNG droplet up to time t. The algorithm to determine 
h{x,t), \x\ < t, for the droplet involves the same permutation as for Hx{t), only the 
partition has to be into increasing subsequences starting from the right. Nevertheless 
h{x, t) equals the longest increasing subsequence (from the left) by the same reasoning 
as in Section ^. As a result we have the alternative expression 

H {t) = max h{t , x) in distribution, (16) 

|x|<i 

where h{x,t) is the height process from the PNG droplet. 
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We note that the maximuin is taken among highly correlated random variables. 
To model their statistics we assume local stationarity and write for the height of the 
droplet h{x,t) = h{0,t) — x^/(t-\/2) + b{x) in the vicinity of a; = 0, where b{x) is a two 
sided random walk with 6(0) = 0, zero mean and the covariance (6(x)^) oc |a:|, as long 
as X = o{t). Therefore the maximum in (^) is achieved inside a region of order t^/^ 
around x = and the modifications between the distributions of H{t) and h{0,t) are 
on the scale t^^^. Our argument suggests that x ¥" ^^ which is indeed supported by 
the numerical simulations of Section ^. In addition it points at a close but non-trivial 
connection between the distributions ^ and x which might pave the way for an analytical 
description of C,. 

4. Universality 

As explained in the Introduction, when starting with a flat substrate of slope u, the 
height at the origin scales as 

H{t) = v{u)t + C{u)t^/^^. (17) 

v{u) and C{u) are non- universal parameters which depend on the slope u and have to 
be computed separately for each model according to the definition in (||) and above. 
However, ^ is universal under the stated initial conditions. 

Clearly it is desirable to have a corresponding scaling form for droplet growth, the 
combining feature being the scale invariance of the macroscopic shape, which means 
that for large / we have 

h{lx,lt) = lh{x,t) (18) 

with h denoting the deterministic profile. For a flat substrate ([T8|) holds trivially, 
whereas for a droplet the shape is determined by the Wulff construction (see for example 
pl|). In other words, the deterministic growth equation 

dMx,t) = v{dMx,t)) (19) 

with the inclination dependent growth velocity v allows only for two types of self-affine 
solutions, the flat surface with slope u, h{x,t) = tv{u) + xu and the droplet 

h{ct,t)=tf{c), (20) 

where / solves the equation /(c) — cf\c) = f (/'(c)). The parameter c is related to the 
local slope u = f'{c) by c = —v'{u). In view of ([17D we are therefore led to the universal 
scaling law for the height distribution of a growing droplet 

Met, t) =t/(c) + sign(A)(i|A|A2t)i/3^. (21) 

Here x is a random variable with TW distribution. 

With the rigorous result (|T3p we are in the position to check the scaling form (pi|) 
for the PNG droplet. In the steady states of the PNG model the left and right island 
edges (steps) have an independent Poisson distribution whose strength is determined 
by the average slope ||22|, ^]. With this information one obtains, in our units, the slope 
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dependent growth velocity v{u) = \J2 + m^, yielding A = 2/(2 + m^)"^/^, and the static 
roughness amplitude A = v{u). The relation c = —u/yl + v? is invertible, therefore one 
can write XA^ = a/2(1 — c^). Comparing with (|T3|) we see that (^Tj) is indeed satisfied 
for the PNG droplet. In addition, with the definition (JT^) of the TW distribution, the 
dimensionless scale factor, denoted by a in Equation (^), is fixed to be |. 

In a remarkable paper |^ Johansson treats the totally asymmetric simple exclusion 
process (TASEP) by a similar but somewhat more involved mapping to increasing 
subsequences in generalized permutations. Interpreted as a one-dimensional discrete 
growth model it describes a crystal occupying initially three quadrants of the plane. 
Further particles, represented as unit squares, may nucleate only if at least two nearest 
neighbors are already occupied. Thus the crystal grows into the forth quadrant starting 



from the origin. It has been known for some time ||25| that on a macroscopic scale the 
shape becomes deterministic. Johansson proves a scaling law similar to (0) with the 
same limiting distribution, but a negative prefactor. In the TASEP particles jump, 
subject to the exclusion constraint, to the right with rate 1. Since the stationary 
measures are Bernoulli, one has v{u) = \{1 — u^), \u\ < 1 and A{u) = \{1 — u^). 



Inserting in (|21|), we recover the result of Johansson, the negative prefactor coming 
from \{u) = —\- This is one of the few instances where universality has been rigorously 
estabhshed on the level of distributions. 

5. Numerics 

In Figure 2 we compare the two distributions of x ^"^^ ^ in a semi- logarithmic plot. The 
TW distribution, which is drawn as a solid line, is obtained by numerically solving (]TB|). 
Its mean is —1.771087 and its variance is 0.813195. For the skewness, {x^)c/{x^)c , 
and the kurtosis, {x^)c/{x^)l, one obtains 0.2241 and 0.09345, respectively. Here (x'^)c 
denotes the k-th cumulant (fully truncated moment). Superimposed are the empirical 
distributions determined at various times, ranging from 1000 to 2000 units, from a 
Monte-Carlo simulation of the PNG droplet. Finite time effects are most pronounced 
for the first moment. At t = 2000 there is still a noticeable deviation of about 3% which 
seems to fade away very slowly. The other quantities agree with the exact values within 
the statistical fiuctuations. 

Since for the distribution of ^ there is no analytic expression available, we performed 
Monte Carlo simulations of the continuum PNG model for fiat and equidistantly stepped 
initial conditions. We used slopes 0, | and 1, up to only moderate times t < 2500 on an 
effectively infinite substrate. The number of independent runs were around 10 000 for 
each slope, but due to translation invariance the number of sample points is about 10^ 
for each recorded time. From our data we estimate the mean, variance, skewness, and 
kurtosis of $, as —0.76(1), 0.637(3), 0.29(1), and 0.16(1), respectively. In view of the 
slow convergence of the first moment for the PNG droplet the estimate for (^) has to 
be taken with care, since it might not yet have reached its true asymptotic value. Our 
values are consistent with earlier studies of C, for various growth models |P, [l^ . 
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Figure 2. Semi-logarithmic plot of the scaled height distributions. The solid line is 
the Tracy- Widom distribution. Superimposed are Monte Carlo results for the PNG 
droplet. The right data points are Monte Carlo results for the flat case as explained 
in the text. 



So far we started the growth process with a deterministic height profile. An obvious 
variant is to start with a random configuration. The most natural choice is to take 
a stationary height statistics obtained in the long time limit. If h{x,t) denotes this 
stationary growth process (i.e. the statistics of dr^h^x, t) depends only on space-time 
differences), then as before 

h{v\u)t, t) - /i(0, 0) ~ v{u)t + C{u)t^'^r] (22) 

for large t. Along a trajectory with slope different from v'{u) one has ordinary diffusive 
scaling as t^^^ due to the static roughness of the height profile. Our simulations (not 
shown) yield, the mean being zero by definition, 1.476(7), 0.375(5), and 0.31(1) for the 
variance, skewness, and kurtosis of r/, which again is consistent with [0] and [HI 



6. Conclusions 



We have identified three distinct universal distributions: growth with noisy stationary 
initial data, and curved and fiat growth from deterministic initial data. In addition we 
have pointed out that the height of the PNG droplet is simply related to the statistics 
of the length of the longest increasing subsequence of a random permutation. This leads 
to a complete characterization of the macroscopic shape and of the fluctuations in the 
height above a single reference point. 

It is somewhat unexpected to have the links one-dimensional surface growth 
to random permutations to random matrices. As a first introduction to the latter 
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we recommend the survey paper by Aldous and Diaconis ||T6[. Longest increasing 



subsequences arise there in the context of sohtaire or patience card games. In patience 
sorting a well shuffled deck of A^ cards, labeled 1, 2, . . . , A^, is dealt into piles according 
to the, admittedly somewhat boring, two rules: (i) A low card is always placed on the 
leftmost pile with a higher card on top. (ii) If (i) cannot be achieved, a new pile is started 
to the right of the previous piles. The connection to the PNG model is immediate. A 
pile of cards corresponds to a height line, the pile size is the number of nucleation events 
on the height line, and, of course, the number of piles is the height. As regards to surface 
growth it would be of great interest to further exploit the random matrix theory in the 
computation of universal properties and to extend this technique to higher dimensions. 
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Addendum 

In the meantime we have noted that the scaling distribution ^ for a flat substrate is the 
limiting distribution of the largest eigenvalue of a GOE distributed random matrix |2^ . 
More precicely, the distribution function of ^ is given by 

F{2^/H) = exp ( - - / (x - t)u{xfdx - - / u{x)dxy (23) 

Its first four cumulants are —0.76007, 0.63805, 0.2935, and 0.1652. Our proof uses the 
symmetrized random permutations studied by Baik and Rains |2^. We are grateful to 
Percy Deift and Craig Tracy for pointing out this reference to us. 
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